Detection and quantification of inverse spin Hall effect from spin pumping in 

permalloy/normal metal bilayers 
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Spin pumping is a mechanism that generates spin currents from ferromagnetic resonance (FMR) 
over macroscopic interfacial areas, thereby enabling sensitive detection of the inverse spin Hall effect 
that transforms spin into charge currents in non-magnetic conductors. Here we study the spin- 
pumping-induced voltages due to the inverse spin Hall effect in permalloy/normal metal bilayers 
integrated into coplanar waveguides for different normal metals and as a function of angle of the 
applied magnetic field direction, as well as microwave frequency and power. We find good agreement 
between experimental data and a theoretical model that includes contributions from anisotropic 
magnetoresistance (AMR) and inverse spin Hall effect (ISHE). The analysis provides consistent 
results over a wide range of experimental conditions as long as the precise magnetization trajectory 
is taken into account. The spin Hall angles for Pt, Pd, Au and Mo were determined with high 
precision to be 0.013 ± 0.002, 0.0064 ± 0.001, 0.0035 ± 0.0003 and -0.0005 ± 0.0001, respectively. 

PACS numbers: 72.25.Rb, 75.47.-m, 76.50.+g 



I. INTRODUCTION 

Information in semiconductor electronic devices and 
data storage technologies is mainly transported and ma- 
nipulated by charge currents. With advancing minia- 
turization, heat dissipation and power consumption be- 
come significant obstacles to further technological ad- 
vances. Alternative technologies that solve or at least 
circumvent these problems are needed. One promising 
candidate to replace existing charge-based technologies 
is based on using spin currents; an effort referred to as 
spintronics. 1 Magnetoelectronic devices employing spin- 
polarized charge currents are already actively in use in 
hard drive read-heads and non-volatile magnetic random 
access memories (MRAMs). 

Pure spin currents that are not accompanied by a 
net charge current, may offer additional advantages in 
applications, 2,3 such as reduced power dissipation, ab- 
sence of stray Oersted fields, and decoupling of spin and 
charge noise. Furthermore, undisturbed by charge trans- 
port, pure spin currents can provide more direct insights 
into the basic physics of spin-dependent effects. Pure 
spin currents can be created by, for instance: (i) Non- 
local electrical injection from ferromagnetic contacts in 
multi-terminal structures; (ii) optical injection using cir- 
cularly polarized light; (iii) spin pumping from a pre- 
cessing ferromagnet; and (iv) spin Hall effect. The last 
possibility is particularly interesting, since ferromagnets 
are not involved^— The spin Hall effect is caused by 
the spin-orbit interactions of defect scattering poten- 
tials or the host electronic structure. The efficiency of 
this spin-charge conversion can be quantified by a single 
material-specific parameter, viz. the spin Hall angle 7, 
which is defined as the ratio of the spin Hall and charge 
conductivities^ and can be measured by magnetotrans- 
port measurements Previous experimental studies 



report quite different 7 values for nominally identical 
materials. For example, for Au a giant 7 = 0.113 was 
reported, 1 — while subsequent experiments found values 
that are one or even two orders of magnitude smaller] 13 : 14 
Similarly, for Pt different experiment s 10 ' 14 : 15 resulted in 
7 values that vary between 0.0037 and 0.08. 

Recently we demonstrated a robust technique 14 to 
measure spin Hall angles with high accuracy in arbitrary 
conductors. Our approach is based on the combination 
of spin pumping, which generates pure spin currents, and 
measurements of electric voltages due to the inverse spin 
Hall effect (ISHE)^ Here we present a detailed discus- 
sion of the measurements in ReffLU and examine the va- 
lidity of the theoretical model used to describe the volt- 
ages induced in the Ni 80 Fe2o(Py) /normal metal (N) bi- 
layers. In particular we measure the ISHE voltage as a 
function of angle of the applied magnetic field, and mi- 
crowave frequency and power. We find excellent agree- 
ment between model calculations and experimental re- 
sults. Accounting for the proper magnetization trajec- 
tory is important for a quantitative interpretation of the 
results. Good agreement between the theoretical model 
and experiments for a wide range of controlled experi- 
mental parameters implies that our approach is robust 
and can be used to determine the magnitude and sign of 
spin Hall effects in more conductors than included in the 
present study. 



II. COUPLING BETWEEN SPIN AND CHARGE 
CURRENTS 

A. Spin pumping in Py/N bilayers 

Spin pumping generates pure spin currents in normal 
metals (N), when they are in contact with a ferromag- 
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FIG. 1: (Color online) Schematic model of spin pumping in 
Py/N bilayer. s shows that the polarization of the spin cur- 
rent is oscillating in time with the frequency of the magne- 
tization precession. The polarization of the spin current is 
perpendicular to the instantaneous magnetization direction 
m, and the rate of magnetization change drh/dt. 



net with time-dependent magnetization induced, e.g., by 
ferromagnetic resonance (FMR)»i2r— The instantaneous 
spin-pumping current j° at the Py/N interface is given 
by^ 
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where rh is the unit vector of the magnetization, s is 
the unit vector of the spin current polarization in N, and 
Re(g'^) is the real part of the spin-mixing conductance. 
The spin current generated by spin pumping is polarized 
perpendicular to the instantaneous magnetization direc- 
tion m and its time derivative drh/dt (see Fig. [T]). Note 
that this spin current always has a polarization compo- 
nent along Hdc and propagates into N normal to the in- 
terface. 

The spin current generated at the Py/N interface ac- 
cumulates a spin density inside the N layer. In the 
ballistic limit (i.e., no spin relaxation in N) the spin cur- 
rent reaching the N/vacuum interface is fully reflected 
and reabsorbed upon returning to the Py/N interface, 
without influencing the magnetization dynamics of the 
bilayer system. 

In real systems, pure spin currents are not conserved, 
since spins relax over length scales given by the spin dif- 
fusion length X s d in N, and the accumulated spin density 
moves across the N layer via spin diffusion limited by 
momentum scattering (leading to electrical resistance) 
and spin-flip scattering (leading to loss of spin angular 
momentum) by spin orbit coupling or magnetic impuri- 
ties. The spin diffusion equation describes the dissipative 
propagation of the spin accumulation (difference in local 
electrochemical potentials of up and down spins) /In in 
the N layer: 
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where u> is the angular frequency, r a f is the spin-flip time, 
z is the coordinate normal to the interface, and D = 
VpTei/3 is the electron diffusion constant, with r e ; the 
electron momentum relaxation time^ The solutions of 
Eq. (|2|) depend on the boundary conditions. For a single 
magnetic layer structure Py/N the boundary condition 
at the Py/N interface is given by^ 
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while for the outer interface we use the free magnetic 
moment condition (full spin current reflection) 
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Equations (J2])-(j4|) can be solved analytically to yield the 
decay of the spin accumulation as a function of the dis- 
tance from the Py/N interface. This decay results in 
spin accumulation profile in the N layer, which decays as 
a function of the distance from the interface, thus driving 
a spin current with a dc contribution: 



3s(z)=f s 



nh[(z - t N )/X sd ] 
smh(t N /X sd ) 



(5) 



where tjy is the thickness of the N layer. The spin accu- 
mulation in N gives rise to spin backflow into the fer- 
romagnet, which effectively reduces the spin pumping 
current, which can be accounted for by replacing g'^ in 
Eq. ([I]) with an effective spin mixing conductance g^j^Mi 
In FMR experiments the absorption of the microwave 
field that excites the magnetization is monitored. The 
magnetization dynamics in ferromagnetic films can be 
described by the Landau-Lifshitz-Gilbert (LLG) equation 
of motion: 
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where "f g = ge/2mc is the absolute value of the gyromag- 
netic ratio, and olq is the dimensionless Gilbert damping 
parameter. The first term on the right-hand side repre- 
sents the precessional torque due to the effective internal 
field H e ff, which for the case of permalloy with small 
anisotropy is approximately equal to the externally ap- 
plied magnetic field H dc - The second term in Eq. ^ rep- 
resents the Gilbert damping torque i 25 ' 26 The spin pump- 
ing can be accounted for in the LLG equation of motion 
by adding a spin pumping contribution a sp to ac, i.e., 
the effective damping becomes at e ff = a G + a sp- The 
damping can be quantified by measuring the FMR line 
width AH, half width at half maximum (HWHM), of 
the imaginary part of the rf susceptibility \ i which is 
commonly measured at a constant microwave frequency 
by sweeping the dc magnetic field Hd c - In case of Gilbert 
damping, AH depends linearly on the microwave angu- 
lar frequency uit, i.e., AH = a e ffU>f/~/ g . The difference 
in the damping parameter, determined by the FMR line 



width, for samples without capping layer and samples in 
which the capping N layer is sufficiently thick to fully 
dissipate the pumped magnetic moment, is attributed 
to the loss of spin momentum in Py due to relaxation 
of the spin accumulation in N. This permits the deter- 
mination of the additional interface damping due to spin 
pumping, 27 which in turn fixes the interfacial spin-mixing 
conductance to: 
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where tp y is the Py layer thickness, M s is the Py sat- 
uration magnetization, and (j,b is the Bohr magneton. 
Note that Eq. (J7J is only applicable when the damping 
is governed by the Gilbert phenomenology or AH ccw/, 
i.e., when inhomogeneous linewidth broadening is neg- 
ligible. Otherwise sjjj can still be determined from 
the additional Gilbert-like damping contribution a sp — 
a Py\N ~ a Pyi where the latter two are obtained from the 
linewidth difference that scales linear with frequency, i.e., 
AH = AH ih + a e ffU)f/ n f g , where AH ih is the sample- 
dependent inhomogeneous linewidth, measured as the 
zero-frequency intercept. 

The dc component of the spin current pumped into N 
is polarized parallel to the equilibrium magnetization and 
has previously been detected via a dc voltage normal to 
the Py/N interfaced Under a simple circular precession 
of the Py magnetization the time-averaging of the spin 
current from Eq. (JT|) for small precession cone angles 9 
reads: 
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In thin magnetic films the trajectory of the magnetization 
precession is not circular but elliptic due to the strong 
demagnetizing fields, which force the magnetization into 
the film-plane. The time-dependent cone angle modi- 
fies the dc component of the pumped spin current by an 
ellipticity correction factor P as derived and measured by 
Ando et al.£^ For an in-plane equilibrium magnetization 
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Equation is a non-monotonic function of Uf, and P 
can become slightly larger than 1, but tends towards 1 
for high frequencies, i.e., large applied fields. 



B. Inverse spin Hall effect 




waveguide 



FIG. 2: (Color online) Experimental setup: (a) Optical image 
of the Py/Pt bilayer integrated into the coplanar waveguide, 
(b) Contacts are added at the end of the bilayer to measure 
the voltage along the waveguide direction. 



edges. Spin pumping generates dc and ac components to 
the spin current: j e JJ c and an rf component transverse 
to the equilibrium magnetization direction. In this pa- 
per we address only the ISHE effect generated by the dc 
component j e JJ c - 

The dc ISHE transverse charge current reads: 



? c * H {z)= 1 {2e/n)jfj c [nx®], 



(10) 



where 7 is the spin Hall angle, n is the unit vector nor- 
mal to the interface and (s) is the polarization vector of 
the dc spin current. For jt L the spin-polarization (s) 
is along the equilibrium magnetization direction in Py. 
The dc electric field lies in N in the plane of the films 
and perpendicular to the equilibrium magnetization of 

p 14,16.30 



III. EXPERIMENTAL RESULTS 

Here we elucidate our method to obtain voltage sig- 
nals due to the ISHE in various Py/N combinations un- 
der FMR conditions, thereby determining the spin Hall 
angle 7 with high accuracy. The measured voltage sig- 
nals scale with the sample length and, therefore, can be 
increased readily by making the samples longer. We iden- 
tify two contributions to the dc voltage: one stems from 
the anisotropic magnetoresistance (AMR) and the sec- 
ond from the ISHE, which can be distinguished by their 
symmetries with respect to the field-offset from the reso- 
nance field. Furthermore, we present a theoretical model 
for the spin Hall angle contribution and test its functional 
dependence of several parameters that can be controlled 
experimentally. 



Spin-orbit coupling or magnetic impurities give rise to 
different scattering directions for electrons with opposite 
spin. In their presence, a spin current in N induces a 
transverse Hall voltage. This ISHE transforms spin cur- 
rents into electrical voltage differences over the sample 



A. Experimental set-up 

The Py/N bilayers were integrated into coplanar 
waveguides with additional leads in order to measure 
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FIG. 3: (Color online) Experimental data (symbols) for FMR 
peak positions and FMR linewidths of a Py/Pt bilayer are 
shown as a function of rf frequency in (a) and (b), respectively. 
The solid line in (a) represents the fit to Eq. (|lip and results 
in M 3 = 852 G. The solid line in (b) represents a linear fit to 
the linewidth vs. frequency dependence. 

the dc voltage along the sample. This is shown in 
Fig. [5] for a Py/Pt bilayer, with lateral dimensions of 
2.92 mm x 20 /mi and 15-nm thick individual layers. The 
bilayer was prepared by optical lithography, sputter de- 
position, and lift-off on a GaAs substrate. Subsequently, 
we prepared Ag contacts for the voltage measurements, 
covered the whole structure with 100 nm of MgO (for dc 
insulation between bilayer and waveguide), and defined a 
30- /im wide and 200-nm thick Au coplanar waveguide on 
top of the bilayer. Similar samples were prepared with 
Pd, Au and Mo layers replacing Pt. 

The high bandwidth of the coplanar waveguide setup 
enabled us to carry out measurements with microwave 
excitations in the frequency interval of 4-11 GHz. The 
power of the rf excitation was varied from 15 to 150 mW. 
For a given frequency, experiments were carried out as 
a function of external magnetic field Hd c , with an in- 
plane orientation that could be rotated to arbitrary an- 
gles a with respect to the central axis of the coplanar 
wave guide. While the FMR signal was determined from 
the impedance of the waveguide^ the dc voltage was 
measured simultaneously with a lock-in modulation tech- 
nique as a function of H c i c - 



B. FMR measurements 

The FMR frequency vs. peak position for the Py/Pt 
sample is shown in Fig. Fitting the data to the 

Kittel formula: 

{Lo } h g f = H dc (H dc + 4ttM s ), (11) 

results in the saturation magnetization for Py of M s — 
852 G. Figure Elb) shows the FMR line width as a func- 
tion of frequency. The linear behavior of the FMR line 
width indicates that damping in Py is governed by the 
intrinsic Gilbert phenomenology and any extrinsic effects 
are small. 

Figure @] shows FMR spectra for a Py/Pt bilayer and a 
Py single layer at 4- GHz excitation frequency. The FMR 
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FIG. 4: (Color online) Derivative FMR spectra at 4 GHz for 
(a) Py and (b) P/Pt. Solid lines are Lorentzian line shape 
fits. 



peak positions for the two samples are similar. The main 
difference between the spectra is their FMR linewidth. 
The FMR line widths (HWHM) extracted from fits to 
Lorentzian absorption functions are AH Pt / Py — 16.9 Oe 
for Py/Pt and AH Py = 12.9 Oe for Py. The differ- 
ence in AH can be attributed to the loss of pumped 
spin momentum in the Pt layer. The thickness of the Pt 
layer is 15 nm, which is larger than — 10 ± 2 nm< 32 
Thus, all pumped spin momentum is dissipated in the 
Pt layer and we can extract the value of the spin mix- 
ing conductance g^jj: from the increased linewidth. Us- 
ing Eq. ([7]) we calculate a spin mixing conductance 
glj = 2.1 x 10 19 m~ 2 at the Py/Pt interface. This 
experimental value is somewhat smaller than the previ- 
ously reported 2.58 x 10 19 m" 2 . 24 * 33 Cao et alM showed 
that for high power rf excitation, the spin mixing con- 
ductance can be reduced due to the loss of coherent spin 
precession in the fcrromagnet. This could be the case 
here, since the cone angle for the FMR at 4 GHz is rela- 
tively large, and a slightly larger mixing conductance for 
the smaller precession angles at 11 GHz would lead to 
more consistent frequency dependent values of the spin 
Hall angles as discussed below. 



C. dc Voltage due to ISHE and AMR effects 

Figure [5] shows the dc voltage measured along the sam- 
ples with an external field applied at 45 degrees from 
the coplanar waveguide axis. For the Py/Pt sample we 
observe a resonant increase in the dc voltage along the 
sample at the FMR position. However, the lineshape is 
complex: below the resonance field the voltage is nega- 
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FIG. 5: Voltage measured along the samples vs. field Hdc for 
Py and Py/Pt at 4 GHz is shown with symbols in (a) and (b) 
respectively. Only the AMR contribution is present in the 
Py sample. Solid line in (a) shows a fit to Eq. (|12[) . Both 
AMR and ISHE effects are observed in the Py/Pt. Dotted 
and dash-dotted lines show the AMR and ISHE contribu- 
tions, respectively, which are extracted from fitting the data 
to Eqs. (|12p and (|18[) : the solid line in (b) shows the combined 
fit for the Py/Pt sample. 



tive, it changes sign just below the FMR resonance field, 
and has a positive tail in the high field region. In con- 
trast, the single layer Py sample, which is not affected by 
spin pumping, shows a voltage signal that is purely an- 
tisymmetric with respect to the FMR position and thus 
mirrors the derivative FMR signal shown in Fig. |4jb). 
The voltage generated by the ISHE depends only on the 
cone angle 9 of the magnetization precession [see Eq. ([8])] 
and thus must be symmetric with respect to the FMR 
resonance position. This means that the voltage mea- 
sured in the Py/Pt sample has two contributions: (i) a 
symmetric signal due to the ISHE, and (ii) an antisym- 
metric signal of the same origin as that in the Py control 
sample. 

The antisymmetric voltages observed in both Py and 
Py/Pt depend on the cone angle 9 of the magnetization 
dynamics, since they vary rapidly around the FMR res- 



onance position. This suggests that the antisymmetric 
signal originates from the AMR.— ~— Although the MgO 
provides dc insulation between the sample and the waveg- 
uide, see Fig. Efa), the strong capacitive coupling allows 
some leakage of the rf driving current I r f = I^jsinojft 
into the sample, I r f,s, which flows along the waveguide 
direction. Its magnitude can be estimated from the ratio 
between the waveguide resistance R wg and the sample 
resistance Rs- Irf,s = IrfRwg/ Rs, since the capacitive 
coupling impedance is negligible for ujf > 4 GHz. Fur- 
thermore, due to the strong capacitive coupling (50 pF) 
between sample and wave guide, both rf currents in the 
sample and waveguide are for all practical purposes in 
phase, i.e., the relative phase shift is expected to be at 
most 10~ 3 7r. Indeed, experiments with single layer Py 
samples^ and with a MgO layer inserted between the 
Py and Pt layers^ are consistent with a pure AMR sig- 
nal, as described below, without any appreciable phase 
shift. 

The precessing magnetization in the Py [see Fig. [SJb)] 
results in a time-dependent Rs[ip(t)] = Ro — 
ARamr sin 2 ?p(t) due to the AMR given by ARamr, 
where Rq is the sample resistance with the magnetization 
along the waveguide axis, and ijj is the angle between the 
instantaneous magnetization fh and the waveguide axis 
[see Fig. EJb)].— ARamr can be experimentally deter- 
mined by static magnetoresistance measurements under 
rotation of an in-plane field sufficiently large to saturate 
the magnetization. Since the AMR contribution to the 
resistance oscillates at the same frequency as the rf cur- 
rent, but phase-shifted, a homodyne dc voltage develops 
and is given by>H 
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Costco , (12) 



where tpo is the phase angle between magnetization pre- 
cession and driving rf field, and the relation between 9, a 
and tp is illustrated in Fig.^b). Well below the FMR res- 
onance the phase angle <po is zero, it becomes tt/2 at the 
peak, and is it far above the resonance.— Thus, cos <po 
changes sign upon going through the resonance, which 
gives rise to an antisymmetric Vamr, as is observed in 
both the Py and Py/Pt samples. Following Guan et al£2- 
we calculate the cone angle 9 and sin ipo as a function of 
the applied field Hdc, FMR resonance field H r , FMR 
linewidth AH and rf driving field h r f. 
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The anisotropic magnetoresistance was determined by 
dc magnetoresistance measurements with fields applied 
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FIG. 6: (Color online) Splitting of rf current due to capaci- 
tive coupling is schematically shown in (a) together with the 
directions of the applied dc magnetic field Hdc and the rf 
driving field h r f with respect to the bilayer and waveguide, 
(b) Schematic of rh precessing in Py. m precesses around its 
equilibrium direction given by Hdc at the driving frequency ojf 
and with a phase delay ipo with respect to h r f. a is the angle 
between Hdc and the waveguide axis (along y), 6 is the cone 
angle described by rh and ip is the angle between fh and the 
waveguide axis. Due to the strong capacitive coupling part of 
I r f flows through the Py given by I r f.s- (c) Geometry of the 
dc component of the pumped spin current with polarization 
direction (s) along the equilibrium magnetization direction 
win. The charge current due to ISHE j ! c SH is orthogonal to 
the spin current direction (normal to the interface) and (s). 
The voltage due to the ISHE is measured along y (waveguide 
axis). Solid arrows indicate the spin accumulation inside N, 
which decays with the spin diffusion length A st j. 



along the hard-axis as ARamb. = 0.95%. This allows us 
to fit the Py data [see Fig.[5fa)] with only one adjustable 
parameter h r f — 4.5 Oe using Eqs. (| 121 TT4|) . 

In order to understand the symmetric contribution to 
the Py/Pt voltage data we have to include an additional 
voltage due to the ISHE. In principle an inductive cou- 
pling (if any) could result in a symmetric voltage con- 
tribution to the signal. However, this type of coupling 
is unlikely in our samples due to the fact that sample 
and transmission line are prepared as a stack with a thin 
insulator in the middle. Furthermore, our recent work 39 
showed that if spin pumping is suppressed by inserting a 



3-run MgO layer at the Py/Pt interface, then the sym- 
metric part of the voltage vanishes. This unambiguously 
shows that the symmetric part of the measured voltage 
is related to spin accumulation in N, which appears due 
to the ISHE. The absence of a symmetric contribution 
for Py alone also suggests that inductive effects are neg- 
ligible. 

In an open circuit, an electric field E is generated lead- 
ing to a total current density 



j c (z) = j c (z) (x cos a + y sin a) + <jE, 



(15) 



where a is the charge conductivity and x, y are defined 
in Fig. [6] Since there is no current flowing in the open- 
circuit 



In 



j c (z)dz = 0. 



(16) 



When the wire is much longer than thick, the electric 
field in the wire is constant. On the other hand, voltage 
generation occurs only in the Pt layer (more precisely in 
a skin depth of the spin- flip length in which the ISHE emf 
is generated), while the Py layer acts as a short, which 
decreases the voltage difference at the sample terminals. 
Solving the system of Eqs. (|T5"|) and (TT6")) . we obtain the 
component of the electric field along the measurement 
direction y as 
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Pgljf sin a sm2 #7 ew f^sd 
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where ctjv and ap y are the charge conductivities in the N 
layer (e.g., Pt) and Py, and tjv and tp y are the thicknesses 
of the N and Py layers. Using Eq. (flTf we calculate the 
voltage due to the ISHE generated along the sample with 
length L: 
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(18) 

Note that this voltage is proportional to L and, thus, 
sufficiently large voltage signals can be measured even 
for small 7 values by increasing the sample length. Fur- 
thermore, note that for the case of the normal layer thick- 
ness 1 7v being comparable to the spin diffusion length X s d 
the measured voltage depends only very weakly on either 
value, since iAr/A sc itanh(ijv/2A S£ i) is approximately con- 
stant. 

One of the input parameters in Vish is the ellipticity 
correction factor P. At 4 GHz excitation, FMR occurs 
at Hdc ~ 200 Oe. Therefore, the magnetization preces- 
sion trajectory is highly elliptical and a correction to the 
dc voltage component due to the ellipticity is significant. 
Figure [7] shows P as a function of microwave frequency as 
calculated using Eq. ©. In the range from 4 to 13 GHz, 
P changes almost by a factor of 3 and, therefore, has 
to be taken into account. At frequencies above 10 GHz 
P becomes larger than 1 and reaches a maximum value 
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FIG. 7: Elliptical precession trajectory results in a time- 
dependent cone angle of magnetization precession, that modi- 
fies the dc component of pumped spin current. The ellipticity 
correction factor P for the dc component of spin current is 
calculated as a function of microwave frequency according to 
Eq. ©. 



of 1.3 at w 28 GHz before it slowly decreases towards 1 
for higher frequencies. This means that the most effec- 
tive pumping of dc component of spin current is achieved 
not for circular precession, but rather for some elliptical 
trajectory of magnetization precession. 

We used Eqs. (fT5|) and (IT2"j) to fit the voltage measured 
for the Py/Pt sample [see the solid line in Fig.JSfb)]. The 
dashed and dotted lines in Fig. EJb) are the AMR and 
ISHE contributions, respectively. By using a literature 
value for Pt of X s d = 10 ± 2 nm^ the only remaining 
adjustable parameters are the rf driving field h r f and 
the spin Hall angle 7 « 0.011 ± 0.002. Note that through 
the cone angle 9, h r f affects both the AMR and ISHE 
contributions. In fact, as seen from the fit to the control 
Py sample, h r f is already determined by the negative and 
positive tails of the AMR part. 

We carried out additional measurements of the spin 
Hall angle as a function of the microwave frequency. 
Since the spin pumping is proportional to the time deriva- 
tive of the magnetization [as manifested by the factor w f 
in Eq. (jTSJ] and the ellipticity correction factor P, which 
increases with frequency (see Fig.[7J), the voltages due to 
the ISHE are expected to increase at higher microwave 
frequencies. However, the cone angle of magnetization 
precession for a constant power of rf excitation decreases 
due to higher resonance fields. Since spin pumping is pro- 
portional to sin 2 an overall decrease in the voltage due 
to the ISHE is observed. However, the relative strength of 
the antisymmetric (AMR) and symmetric (ISHE) parts 
of the signal change since the AMR decreases faster than 
the ISHE contribution as a function of the frequency. 
This effect is illustrated in Figs.[S](a) and (b). FigureJS^c) 
shows the values of the spin Hall angle 7 extracted from 
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FIG. 8: (Color online) (a) and (b) show voltages measured 
at 4 and 11 GHz. Voltages measured at 11 GHz are smaller 
due to a decreased precession cone angle. Note that the ratio 
between the ISHE and AMR contributions change due to a 
faster decrease of the AMR voltage at high frequencies, (c) 
Spin Hall angle in Pt as a function of frequency. The slightly 
decreased values at lower frequencies may be due to non-linear 
effects and a concomitant decrease of spin pumping at large 
angles of magnetization precession. 



the fits, which are essentially constant for all frequen- 
cies, except for a slight decrease of 7 at lower frequen- 
cies. Equation © is strictly valid only for small preces- 
sion cones and constant power of rf excitation. From 
the fitting of the data we can extract the cone angles of 
the magnetization precession. At 4 GHz the fitted value 
9 w 10° at the resonance, while at 11 GHz w 2.5°. For 
4-GHz excitation non-linear effects may start to play a 
role, possibly slightly changing the estimated value of 7. 
Our model was further tested by varying the angle a of 
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FIG. 9: (Color online) (a) Voltage measured at 4 GHz as a 
function of angle a of the external magnetic field with respect 
to the coplanar waveguide axis. Experimental data and fits 
are shown with symbols and solid lines, respectively, (b) Spin 
Hall angle extracted from the fits. The theoretical model 
correctly takes into account the angular dependence for the 
ISHE and AMR contributions. 
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FIG. 10: (Color online) Power dependence of the symmetric 
ISHE voltage contribution measured at 4 GHz. The inset 
shows that the maximum of the ISHE signal is linear vs. rf 
excitation power. The highest power (150 mW) deviates from 
the linear behavior due to excitation of non-uniform modes. 
This deviation is also observed in the corresponding FMR 
spectra. 



TABLE I: Spin Hall angle 7 determined using X 3 d and ctjv 



IT 

Normal metal 


)im naxa 1 
X 3d (nm) 


leasuren nr. i i i,n? 
ajv l/(yim) 


7 


Pt 




10±2 


(2.4±0.2)xl0° 


0.013±0.002 


Pd 




15±4 


(4.0±0.2)xl0 6 


0.0064±0.001 


Au 




35±3 


(2.52±0.13)xl0 7 


0.0035±0.0003 


Mo 




35±3 


(4.66±0.23)xl0 6 


-0.0005±0.0001 



the applied field with respect to the microwave transmis- 
sion line. Note that both Eqs. (fT2|) and (fT8|) have besides 
the explicit dependence on a an additional dependence 
through the implicit cv-dependence of 9 given by Eq. (|13[) . 
For small cone angles 9 this results in both Vamr and 
Vish being proportional to sin a cos 2 a. The dependence 
on the dc magnetic field direction is shown in Fig. |H1 The 
measured voltage profile is consistent with the theoretical 
model, and results in a consistently constant fitted value 
of the spin Hall angle in Pt. Due to the specific geome- 
try of the sample we were not able to measure at angles 
close to a — 90°, at which the magnetization dynamics 
cannot be excited, because the component of h r f perpen- 
dicular to the magnetization vanishes and FMR cannot 
be excited. But in the range of angles from —5° to 45°, 
excellent agreement between experiment and theory was 
achieved. 

The other adjustable parameter in our measurements 



is the power of the microwave excitation. The rf mi- 
crowave field amplitude increases as a square root of the 
power. According to Eq. fll~3|) the cone angle 9 of the 
magnetization precession increases linearly with driving 
field. The voltage due to the ISHE is quadratic in 9 and, 
thus, is expected to be proportional to the power. After 
fitting the data, we extracted the symmetric part due to 
the ISHE, which is shown in Fig. [TO] The maximum val- 
ues of the measured voltage depend linearly on power, 
as expected by theory, except for the highest power of 
about 150 mW, at which the system is driven into the 
non-linear regime. We observe a deviation of the FMR 
peak position as a well as a deviation of the FMR spec- 
tra from a Lorentzian shape. It is known that at high 
rf power other modes beside the uniform FMR mode 
are excited. In this case one expects a substantial line 
broadening and even saturation of the FMR absorption, 
as observed. 
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FIG. 11: (Color online) Voltages measured at 11 GHz for (a) 
Py/Pd, (b) Py/Au, and (c) Py/Mo. Shown are data (sym- 
bols), combined fits (black lines) and individual AMR and 
ISHE contributions, with dotted green and dash-dotted blue 
lines, respectively. Note the opposite sign of the ISHE contri- 
bution for Mo compared to those for Au and Pd. 



spin Hall contribution has the opposite sign. Fitting of 
the data enabled us to extract the values of 7 for Pd, 
Au and Mo (see Table U). The effective mixing conduc- 
tance at the intermetallic Py/N interface is governed by 
N and we adopt the value obtained by the experimen- 
tally measured increased damping in Py/N. Note that 
the determination of 7 furthermore requires <tn and X s d 
as input parameters. <7n was obtained using four-probe 
measurements for all samples. Reported values for A S( j 
vary considerably. We choose literature values for Pt 
and Pd from Refl32l and Au from Reffl9l. and for Mo we 
assumed that X s d is comparable to that for Au. Even 
though this latter assumption may not necessarily hold, 
the sign change is consistent with earlier measurements* 12 
We furthermore note that the 7 values in Table U differ 
from the previously reported ones in Ref(U, where we as- 
sumed circular precession and therefore underestimated 
7 by a factor of roughly 2. In addition, Table fl] is based 
on 11 GHz data, which due to the smaller cone angles is 
less susceptible to deviations stemming from non-linear 
effects, and therefore should be more reliable. 

Our values for 7 are in good agreement with values re- 
ported by Otani et aZi 12 i 40 from measurements in lateral 
spin valves, but conflict with values reported by other 
groups ! 11 ' 15 We note that in lateral spin valves it is im- 
portant to also understand the charge current distribu- 
tion in order to rule out or correct for additional non-local 
voltage contributions, 13 A distinct advantage of our ap- 
proach is that the measured voltage signal scales with 
the sample dimension, and no additional charge current 
is directly applied to the sample that could result in un- 
wanted spurious voltage signals. 

We also gain insights into spin orbit coupling in non- 
magnetic metals, which ultimately give rise to spin Hall 
effects. Even for non-magnetic materials, that are next 
to each other in the periodic table (Pt and Au) the spin 
Hall angle differs almost by a factor of 4. On the other 
hand, Mo has a significantly smaller spin Hall angle with 
opposite sign. This sign change can be rationalized by a 
simple s-d hybridization model and is supported by first- 
principles calculations^ indicating that the spin Hall an- 
gle should be negative for less than half-filled d-bands, 
and positive for more than half-filled ones, consistent 
with our experimental results. Pd in spite of being a 
lighter element than Au has a spin Hall angle which is 
almost 2 times larger. First-principles calculations are 
again consistent with the experimental observation of 7 
being larger for Pt and Pd compared to Au. 42 



D. Spin Hall angle in Pd, Au and Mo 

Since the sample structure in our experiments is just a 
bilayer of Py with the non-magnetic material of interest, 
this technique can be readily applied to determine 7 in 
any conducting material. In Fig. [11] we show voltages 
measured for Py/Pd, Py/Au, and Py/Mo measured at 
11 GHz excitation. The spin Hall contributions in Au 
and Mo are smaller than in Pt, and note that for Mo the 



IV. CONCLUSIONS 

We presented a spin pumping technique that enables 
measuring spin Hall angles in various materials, which 
has clear advantages over standard dc electrical spin in- 
jection in Hall bar microstructures. Our results for Pt, 
Pd, Au and Mo show that spin Hall angles are rather 
small, with the largest value found to be 0.013 for Pt. 
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Our approach provides a uniform spin current across a 
macroscopic sample. The voltage signal from the inverse 
spin Hall effect can readily be increased via the use of 
longer samples, since Vish oc L. Furthermore, by us- 
ing an integrated coplanar waveguide architecture we can 
control parameters, such as the rf driving field distribu- 
tion, microwave frequency and power of rf excitation. 
This enabled a quantitative analysis of the data and a 
test of the theoretical model under various experimental 
conditions. Our model accounts for both the anisotropic 
magnetoresistance and the inverse spin Hall effect con- 
tributions and agrees with experiments for a wide range 
of controllable parameters. We demonstrated the exis- 
tence of symmetric (ISHE) and antisymmetric (AMR) 
voltages and could model the frequency, magnetic field 
direction, and excitation power dependence well. The 
AMR voltage in our experiments originates from capaci- 
tive coupling between the waveguide and the sample and 
is consistent with the parameters characteristic for the 
ferromagnetic resonance. Our method will enable addi- 



tional studies of spin Hall effects in other materials, and, 
therefore, will be useful to further understand the spin- 
orbit coupling mechanism in metals. This is necessary in 
order to develop and optimize the spin Hall effect as a 
method to generate and detect spin currents in various 
circumstances, such as in the spin Seebeck effect^ 
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